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Abstract 

The energy spectrum for high energy 7-rays (EL > 10 MeV) from 
the process pp — ► 77X emitted at 90° in the laboratory frame has been 
measured at an energy below the pion production threshold, namely, 
at 216 MeV. The resulting photon energy spectrum extracted from 
7 — 7 coincidence events consists of a narrow peak at a photon energy 
of about 24 MeV and a relatively broad peak in the energy range of (50 
- 70) MeV. The statistical significances for the narrow and broad peaks 
are 5. 3d and 3.5a, respectively. This behavior of the photon energy 
spectrum is interpreted as a signature of the exotic dibaryon resonance 
d\ with a mass of about 1956 MeV which is assumed to be formed in 
the radiative process pp — > 7 d\ followed by its electromagnetic decay 
via the d\ — > pp'j mode. The experimental spectrum is compared with 
those obtained by means of Monte Carlo simulations. 

1 Introduction 



Direct-production of two hard photons in nucleon-nucleon collisions (./V./V77) 
at intermediate energies, unlike production of single photons (NN bremsstrahlung) , 
is one of those fundamental processes which are still poorly explored both 
theoretically and experimentally. However, a study of this process can yield 
new important information on the underlying mechanisms of the NN in- 
teraction that is complementary to that obtained from investigations of 
other processes. In particular, the NNjj process can be used as a sen- 
sitive probe for experimental verification of the possible existence of NN- 
decoupled nonstrange dibaryon resonances. These are two-baryon states 2 B 
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with zero strangeness and exotic quantum numbers, for which the strong 
decay 2 B — > iViV is either strictly forbidden by the Pauli principle (for the 
states with the isospin / = 1(0) and with an even sum L + S + I, where 
L is the orbital momentum and S is the spin), or is strongly suppressed by 
the isospin selection rules (for the states with I = 2)|], §. Such dibaryon 
states cannot be simple bound systems of two color-singlet nucleons, and a 
proof of their existence would have consequences of fundamental significance 
for the theory of strong interactions. The iViV-decoupled dibaryon states 
are predicted in a series of QCD-inspired mo dels 0, f|, [|. Among the pre- 
dicted dibaryons, there are those with masses both below and above the pion 
production (wNN) threshold. The iViV-decoupled dibaryons with masses be- 
low the tcNN threshold are of special interest, since they may decay mainly 
into the NN'-f state, and consequently, their widths should be very narrow 
(< IkeV). Narrow dibaryon states have been searched for in a number of 
experiments^], but none has provided any convincing evidence for their exis- 
tence. Most of the past searches, however, were limited to dibaryons coupled 
to the NN channel, and, to our knowledge, only a few were dedicated to the 
iViV-decoupled dibaryons|7|, §, [J. At the same time, since the dibaryons we 
are interested in do not couple to the iViV-channel, then, depending on their 
production and decay modes, they may have escaped detection up to now 
and may naturally appear in dedicated experiments. 

If the iViV-decoupled dibaryons exist in nature, then the NN'yy process may 
proceed, at least partly, through the mechanism that directly involves the 
radiative excitation NN — > 7 2 B and decay 2 B — > 'yNN modes of these 
states. In NN collisions at energies below the nNN threshold, these pro- 
duction and decay modes of the iViV-decoupled dibaryon resonances with 
masses Mr < 2m jy + would be unique or dominant. The simplest and 
clear way of revealing them is to measure the photon energy spectrum of 
the reaction NN'yy. The presence of an iViV-decoupled dibaryon resonance 
would reveal itself in this energy spectrum as a narrow line associated with 
the formation of this resonance and a relatively broad peak originating from 
its three-particle decay. In the center-of-mass system, the position of the nar- 
row lme(En) is determined by the energy of colliding nucleons (W = v^i) and 
the mass of this dibaryon resonance as Er = (W 2 — Mj { )/2W. An essential 
feature of the two-photon production in NN collisions at an energy below 
the pion production threshold is that, apart from the resonant mechanism in 
question, there should only be one more source of photon pairs. This is the 
double iViV-bremsstrahlung reaction. But this reaction is expected to play a 
minor role. Indeed, it involves two electromagnetic vertices, so that one may 
expect that the NN / yy-to-NN / y cross section ratio should be of the order of 
the fine structure constant a. However, the cross section for NN'j is already 
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small (for example, the total pp-bremsstrahlung cross section at energies of 
interest is a few fib). 

The process pp — > pp~f~f at an energy below the pion production threshold 
provides an unique possibility of searching for the iViV-decoupled dibaryon 
resonances in the mass region Mr < 2m p + m v with quantum numbers 
I(J P ) = 3 + , etc.) or those with I = 2 (J is the total spin, and P is the 
parity of a dibaryon). The preliminary experimental studies of the reaction 
pp — ► 77 V at an incident proton energy of about 200 MeVpD|, [TT|] showed 
that the photon energy spectrum of this reaction had a peculiar structure 
ranging from about 20 MeV to about 60 MeV. This structure was interpreted 
as an indication of the possible existence of an iVTV-decoupled dibaryon res- 
onance(later called off) that is produced in the process pp — > 7 d* and subse- 
quently decays via the d* — > ppy channel. Unfortunately, a relatively coarse 
energy resolution and low statistics did not allow us to distinguish the nar- 
row 7-peak associated with the d* production from the broad 7-peak due 
to the decay of this resonance and, hence, to determine the resonance mass 
exactly. To get a rough estimate of the d* mass, we admitted that the ex- 
pected narrow 7-peak is positioned at the center of the observed structure 
and thus obtained Mr ~ 1920 MeV. However, if the uncertainty in the posi- 
tion of the narrow 7-peak in the structure is taken into account, the actual 
d\ mass might be considerably different from such a crude estimate. The 
possible existence of this dibaryon resonance was also probed in the proton- 
proton bremsstrahlung data taken by the WASA/Promice Collaboration at 
the CELSIUS accelerator [I2~] for proton energies of 200 and 310 MeV. How- 
ever, the trigger used in those experimental studies to select events was de- 
signed and optimized for investigating the usual pp — > ppy bremsstrahlung. 
The selected events were those corresponding to two simultaneously detected 
protons emerging in the forward direction at angles between 4° and 20° with 
respect to the beam axis, each of which had a kinetic energy exceeding the 
present threshold (40 MeV). The analysis of those data resulted only in up- 
per limits on the dibaryon production cross section in the mass range from 
1900 to 1960 MeV. Therefore, to clarify the situation with the dibaryon reso- 
nance d*, we have decided to measure the energy spectrum of the pp — > ppjj 
reaction more carefully. 



2 Experimental details and event selection 

The experiment was performed using the variable energy proton beam from 
the phasotron at the Joint Institute for Nuclear Research (JINR). The schematic 
layout of the experimental setup is shown in Fig. 1. The pulsed proton beam 
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Figure 1: Experimental layout. Proton beam approaches from the left. M 
is the monitor of the incident beam, T is the target, D\ and D 2 are the 
7-ray detectors, V\ and V 2 are the plastic scintillators, LC\ and LC 2 are the 
lead collimators, PB 1 and PB 2 are the polyethylene bars. At the top right 
of the figure the cross section of the 7-ray detector D\ is shown. C is the 
central detector module and R±, R 2 , .., Rq (ring) are the detector modules 
surrounding the detector module C. 

(bursts of about 7 ns FWHM separated by 70 ns) with an energy of about 
216 MeV and an energy spread of about 1.5% bombarded a liquid hydro- 
gen target (T). The average beam intensity was on the order of 3.6 • 10 8 
protons/s and was monitored by an ionization chamber (M). The liquid 
hydrogen target was a cylindrical cryogenic vessel with a length of about 5 
cm (~ 0.35 g/cm 2 ) and a diameter of 4.5 cm, which had thin entrance and 
exit windows made from 100 fim thick mylar foil. Both 7-quanta of the re- 
action pp r yy were detected by two 7-ray detectors (D x and D 2 ) placed in a 
horizontal plane, symmetrically on either side of the beam at a laboratory 
angle of 90° with respect to the beam direction. The 7-detector D\ was an 
array of seven individual detector modules (see Fig. 1): the central detector 
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module (C) and six detector modules (R±, R 2 , .., Re) surrounding the central 
detector module. Each detector module was a 15 cm thick CsI(Tl) crys- 
tal having a hexagonal cross section with an outer diameter of 10 cm. The 
detector modules Ri, R 2 , .., Re (ring) were designed to measure the energy 
of secondary particles and photons leaving the central detector module C . 
Thus, the energy of each photon detected by the detector D 1 was determined 
by summing up the energy deposited in the central detector module and the 
energies deposited in the detector modules of the ring. Besides, the ring was 
used as an active shielding from cosmic ray muons. The shielding was based 
on the fact that a cosmic ray muon that passed through the detector D\ 
should deposit an energy at least greater than 50 MeV in the ring . The 
energy resolution of a single CsI(Tl) crystal was measured to be about 13% 
for Ej = 15.1 MeV. The second 7-ray detector D 2 was a cylindrical NaI(Tl) 
crystal 15 cm in diameter and 10 cm in length. Compared to the detector 
D 1 , it had a poorer energy resolution for 7-rays with energies E 1 > 10 MeV. 
The detector D 2 was mainly designed to select 7 — 7 coincidence events. To 
reject events induced by charged particles, plastic scintillators (Vi and V 2 ) 
were put in front of each 7-detector. Both 7-ray detectors were placed in 
lead houses with walls 10 cm thick which were in turn surrounded by a 5 cm 
borated parafin shield. The front wall of every lead house had an opening 
that served as a collimator {LC\ and LC 2 ) limiting the angle acceptance of 
the 7-detector. The solid angles covered by the detectors D 1 and D 2 were 
43 msr and 76 msr, respectively. To attenuate prompt particles coming from 
the target, 15 cm polyethylene bars (PBi and PB 2 ) were placed in between 
the target and each 7-detector. 

The electronics associated with the 7-detectors and the plastic scintillators 
provided amplitude and time signals suitable for digitizing by an ADC and 
a TDC. The energy threshold for both the detector module C and the 
7-detector D 2 was set at about 7 MeV. The amplitude signals from the 7- 
detectors were sent to the ADCs. The time signals from the 7-detectors 
and the plastic scintillators were sent to the TDCs. In addition, the radio- 
frequency (RF) signals from the accelerator were also sent to the TDC. Each 
of these signals was related with the time of arrival of the proton beam burst 
at the target. Therefore, the time spectrum of the RF-signals provided the 
means to select those events which correlate with proton bursts. The time 
signal from the detector module C started all the TDCs and opened the 
ADC gates. In parallel, the time signals from the detector module C and 
the 7-detector D 2 were sent to a coincidence circuit to select those events 
which had signals from both the 7-detectors within a 250 ns time interval. 
The data acquisition system was triggered by signals from the coincidence cir- 
cuit and the corresponding events were recorded in the event-by-event mode 
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on the hard disk of the computer. A further selection of 7 — 7 coincidence 
candidate events was done during off-line data processing. The time interval 
of 250 ns was chosen to record two different types of events, "prompt" and 
"delayed". A prompt event is that in which both photons are due to the 
same beam burst. Such events may be attributed to either real 7 — 7 events 
or to random coincidences. A delayed event is that in which photons are due 
to two different beam bursts, and, therefore, it must be attributed to a ran- 
dom coincidence. This type of events permits one to derive the component 
of random coincidences mixed with real prompt 7 — 7 events. 

Data were taken for the target filled with liquid hydrogen and for the 
empty one. The off-line analysis of both data sets was done in exactly the 
same way. The selection procedure for events associated with the process 
PP ~^ 77A included the following main operations: 

- Rejection of events induced by charge particles coming from the target. 
For this purpose, all events having appropriate time information from 
any plastic scintillators were rejected. 

- Rejection of events induced by cosmic ray muons. These events were 
rejected by a requirement that in every event the energy detected by 
the ring should not exceed 30 MeV. The number of events induced by 
cosmic ray muons was found to be insignificant. 

- Selection of events which are caused by proton beam bursts. This pro- 
cedure was done with the help of a time spectrum of RF-signals. An 
example of this spectrum for a measurement with the full target is 
shown in Fig. 2. The time spectrum of RF-signals consisted of a peak 
of useful events superimposed on a background of events due to ran- 
dom coincidences. The time region confining the peak of useful events 
was determined, and events outside of this region were considered as a 
background and were rejected. 

- Selection of events associated with 7 — 7 coincidences. To do this, a 
spectrum of time signals from the detector D2 (coincidence time spec- 
trum) was used. A part of this spectrum for a measurement with the 
full target including the peak of prompt events and one of the peaks of 
delayed events is shown in Fig. 3. The coincidence time spectrum con- 
sisted of the peak of prompt events and a few peaks of delayed events. 
The time region confining the prompt peak was determined, and events 
associated with it were considered as 7 — 7 coincidence events and were 
accepted. 
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Figure 2: Time spectrum of RF-signals. The dashed lines indicate the time 
interval confining the peak of useful events. 

3 Results and discussion 

Data for the full and empty target were taken in two successive runs for 
~31 h and ~21 h, respectively. The integrated luminosity of about 8.5 pb~ x 
was accumulated for the measurement with the full target. The resulting 
energy spectra of prompt 7 — 7 events as a function of the photon energy 
measured by the detector D\ for the full and the empty target are presented 
in Fig. 4. Apart from real 7 — 7 events, these spectra also contain the 
contributions due to random coincidences which were determined by using 
delayed events. Such contributions to the full- and empty-target spectra 
were found to be identical within the statistical errors. The contribution of 
the random coincidence background to the empty-target spectrum of prompt 
7—7 events was less than 20%. The behavior of the empty-target spectrum is 
determined by photon pairs produced in collisions of protons with the mylar 
windows of the target vessel. Since carbon is the major component of mylar, 
we also performed the measurement with the carbon target to make sure 
that the observed shape of the empty-target spectrum is due to the process 
12 C(p, 77)X. It was found that the empty- and carbon-target spectra have 
the same shape. 

The photon energy spectrum of the process pp — > 77X obtained after 
subtraction of the empty target contribution is shown in Fig. 5. As can be 
seen, this spectrum consists of a narrow peak at a photon energy of about 
24 MeV and a relatively broad peak in the energy range from about 50 
to about 70 MeV. The statistical significance expressed in the number of 
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Figure 3: Coincidence time spectrum. The dashed lines indicate the time interval 
confining the peak of prompt events. 



standard deviations (N.S.D.) was calculated for both the sharp and the broad 
peak using the expression N.S.D. = Q^iV-^i l°i V "* 2 ) 1 ^ 2 ) where 

N- ff and Oi are the number of events and the error bar at the energy 
respectively. As a result, we found that the number of standard deviations 
is equal to 5.3 for the narrow peak and to 3.5 for the broad peak. The width 
(FWHM) of the narrow peak was found to be about 8 MeV. This width is 
comparable with that of the energy resolution of the experimental setup. It 
is important to note that the number of events in the narrow peak (~ 130) 
is approximately equal to half the total number of events in the spectrum (~ 
250). The observed behavior of the photon energy spectrum agrees with a 
characteristic signature of the sought dibaryon resonance d\ that is formed 
and decays in the radiative process pp — >• 7 d\ — > pp'yy. In that case the 
narrow peak should be attributed to the formation of this dibaryon, while 
the broad peak should be assigned to its three-particle decay Using the 
value for the energy of the narrow peak En ~ 24 MeV, we obtained the 
d* mass Mr ~ 1956 MeV. The uncertainty of the given value is primarily 
due to the uncertainties in the value of the proton beam momentum and the 
energy of the narrow 7-peak. We estimate the total uncertainty in the d* 
mass to be ±6 MeV. In such a way we obtain Mr = 1956 ± 6 MeV. This 
result is inconsistent with the rough estimate of the d\ mass Mr ~ 1920 MeV 



we obtained earlier [|H], 11]. The results of the present experiment enable us 



to estimate the differential cross section of the resonance production of two 
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Figure 4: The photon energy spectra of prompt 7 — 7 events for data taken with 
the full target (filled circles) and the empty target (open circles) . 

photons in the process pp — > 7 d\ — > ppjj by the formula: 

d 2 a N eff 

aa = — (1) 

dn 1( in 2 LAQ 1 An2 , 

where N e ** is the total number of events observed, L is the luminosity, 
and AQ 1 and AQ 2 are the solid angles covered by the 7-ray detectors. We 
found that at an energy of 216 MeV this cross section for photons emitted 
symmetrically at 9i a b = ±90° is ~ 9 nb/sr 2 . 

Having assumed that the pp — > 7 d\ — > pp'yy process with the d\ mass of 
1956 MeV is the only mechanism of the reaction pp — * PJ577, we calculated 
the photon energy spectra of this reaction for a proton energy of 216 MeV. It 
was also assumed that the radiative decay of the d* is a dipole E1(M1) transi- 
tion from the two-baryon resonance state to a pp-state in the continuum. The 
calculations were carried out with the help of Monte Carlo simulations which 
included the geometry and the energy resolution of the actual experimental 
setup. The photon energy spectra were calculated for two different scenarios 
of the d\ decay. The difference between them was that one of these scenarios 
took into account the final state interaction (FSI) of two outgoing protons 
whereas in the other that interaction was switched off. Each of the scenar- 
ios imposed some restrictions on possible quantum numbers of the dibaryon 
state in question. The scenario including the FSI implies that the final pp- 
system is in the singlet 1 S state and consequently it should take place, in 
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Figure 5: Experimentally observed energy spectrum for photons from the pp'yy 
process and energy spectra for photons from the process pp — > jd* — > jjpp calcu- 
lated with the help of Monte Carlo simulations for two d\ decay scenarios: without 
the FSI(solid line) and with the FSI(dashed line). The solid and dashed lines are 
overlapped in the energy ranges E 1 < 30 MeV and E 1 > 80 MeV. 



particular, for the isovector 1 + dibaryon state (the simplest exotic quantum 

numbers), namely, 1 + + . Moreover such a scenario can take place for 
any isotensor dibaryon state with the exception of the + or CF state. At the 
same time, the scenario in which the FSI is switched off, is most likely to 
occur for the isotensor 0^ dibaryon state. The spectra calculated for these 
two decay scenarios and normalized to the total number of 7 — 7 events ob- 
served in the present experiment are shown in Fig. 5. Comparison of these 
spectra with the experimental spectrum indicates that both the calculated 
spectra are in reasonable agreement with the experimental one within ex- 
perimental uncertainties. In other words, the statistics of the experiment is 
insufficient to draw any firm conclusions in favor of one of these scenarios 
and thereby to limit possible quantum numbers of the observed dibaryon 
state. In this connection, we would like to emphasize that the experimental 
data on the process n~d — > 77X can shed some light on the properties of 
the dibaryon state we have observed. The point is that the d\ dibaryon can 
in principle be excited in this process via the mode n~d — > ^d\ — > 777171. 
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However, excitation of the dibaryon d\ with isospin I = 1 is allowed by the 
isospin selection rules, whereas excitation of the isotensor (J = 2) resonance 
is expected to be strongly suppressed |14|, [15| . Therefore, the preliminary 
results of the study of the reaction 7i~d — > 77X carried out by the RMC 
Collaboration at TRIUMF|13| indicating no enhancement of the two-photon 
yield in the process {jt~d) atom — > jjnn compared to the usual one-nucleon 
mechanism ir~p — > 7771, on the proton bound in the deuteron may imply that 
the d\ dibaryon has isospin 2. 

Finally, as we have already mentioned in the Introduction, the presence 
of the pp — ► 7 d*(1920) — ► 77PP process was sought in the Uppsala pp- 
bremsstrahlung data||12| . Based on the detailed data analysis the authors fT2| 
conclude that they have found no evidence for a dibaryon in the mass range 
from 1900 to 1960 MeV. There is, however, one important feature of the pro- 
cess in question which was omitted in Ref. ||12]|. It is related to the angular 
distribution of protons from the dibaryon decay d* — > pp'y. The authors JT2] 
assumed that such a distribution is isotropic in the center-of-mass system 
of the reaction. However, that is not the case. To understand the angular 
distribution of the protons from the process pp — > 7 o?*(1956) — ► jjpp we 
performed phase space calculations using the CERN codejnj. The calcula- 
tions were done for two possible decay modes of the dibaryon. One is the 
dibaryon decay to the singlet 1 S two-proton state ( 1 So{pp}). It is allowed, 
for instance, for d* with quantum numbers 1(1 + ) or 2(1 ± ). The second im- 
plies that the dibaryon decays to any other two-proton state. It was also 
assumed that the 7-radiation from the given process has a dipole character. 
The results of the calculations (in the laboratory) for an incident energy of 
310 MeV and for the case when the polar angle of observed photons was 
constrained to the range between 30° and 90° (7-ray detector arrangement 
in the Uppsala experiment) are shown in Fig. 6. The same calculations done 
at 200 MeV led to similar results. As can be seen from Fig.6, the protons 
from the process pp — > 7 d*(1956) — > 77 1 So{pp} are mainly concentrated in 
a narrow angular cone at 6i a b < 4°. Since in the Uppsala experiment events 
were only accepted if the polar angle of each proton fell within a forward 
cone between 4° and 20°, the dominant portion of events associated with the 
given process was beyond this angular range and, hence, were not present in 
the analyzed data[]T2|. In other words, if we assume that the d*(1956) has 
quantum numbers such as 1(1 + ) or 2(1 ± ), then these data can contain only 
such events from the process pp — > 7 o^(1956) — > 77 g 1 S {pp} in which both 
final protons are associated with the tail of the angular distribution presented 
in Fig.6. Our naive estimations have shown that the fraction of these events 
in the data|12[ for 200 MeV and 310 MeV with respect to the correspond- 
ing total number of events from the process of interest are approximately 
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Figure 6: The angular distribution of the protons (in the laboratory) from the 
process pp — > 7 6^(1956) — > 77PP calculated with FOWL for the case when the 
polar angle of observed photons is constrained to the range between 30° and 90°. 
The solid line corresponds to four-body phase space with 1 Sq FSI and the dashed 
line to that without the FSI between the outgoing protons. 



smaller than 20% and 2%, respectively. These estimates together with the 
differential cross section ( d Q^n 2 ) obtained in the present experiment permit- 
ted us, in turn, to estimate the expected number of events in the Uppsala 
experiment. The corresponding values for the measurements at 200 MeV and 
310 MeV are equal to ~2 and ~130, respectively. It was assumed in those 
calculations that the energy dependence of the cross section is E\ (dipole 
radiation), where E R — (s — M^)/2 v /s. Our estimate for 200 MeV is in 
agreement with the corresponding null result obtained in the Uppsala mea- 
surements. As for the result at 310 MeV, we would like to note the following. 
First, the experimental sensitivity of the data[12] at 310 MeV to the process 
in question was limited by a large background stemming from 7r° decays. 
We found that cuts applied in Ref.|jl2f to remove this background reduced 
the number of dibaryon candidates observed by a factor of ~3. Second, the 
presence of a narrow dibaryon should appear in the ppy invariant mass spec- 
trum as a sharp peak and a wider one (mirror) that arises because of the 
PP7 invariant mass calculation with a photon originating from a dibaryon 
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production. The authors [15] believed that inside a range of dibaryon masses 
from 1900 to 1960 MeV these peaks would be well separated from one an- 
other. Therefore, cLS cL SI gnal for dibaryon production they hoped to observe 
a sharp peak ~5 MeV wide. To find the invariant mass distribution of pp^ 
events from the pp — > 7 d*(1956) — > 'yypp process at 310 MeV we performed 
a Monte Carlo simulation of the WASA/PROMICE experimental setup. It 
appeared that a broad mirror and a sharp peak merge in this distribution and 
it consists of one wide peak extending from ~1943 to ~1968 MeV. We think 
that in Ref.|l2| such a peak could be overlooked. Of course, the estimates 
obtained for the expected number of events in the Uppsala experiment are 
very approximate. To make more precise estimates it is desirable to have a 
consistent description of the pp — > jd* — > jjpp process. 



4 Conclusion 

The 7-ray energy spectrum for the pp —>■ 77V reaction at a proton energy 
below the pion production threshold has been measured for the first time. 
The spectrum measured at an energy of about 216 MeV for coincident pho- 
tons emitted at an angle of 90° in the laboratory frame clearly evidences the 
existence of the AA-decoupled dibaryon resonance d\ with a mass of 1956±6 
MeV that is formed and decays in the process pp — > 7 d\ — > ppjj- The data 
we have obtained, however, are still incomplete, and additional careful stud- 
ies of the reaction pp — > ppj^y are needed to get proper parameters (mass, 
width, spin, etc.) of the observed dibaryon state. 
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